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VI. Semiclassical laser Theory 23 

Laser Self-Consistency Equations 

Lamb's theory of laser operation provides a very powerful means for interpreting and 
predicting complex time dependent behavior without invoking all the intricacy of a full 
quantum mechanical theory. It is semiclassical, self-consistent theory in the following 
sense: 



1 



quantum 
theory 



statistic 
summation 



macroscopic 
Maxwell theory 



self-consistency 



Suppose that the field in the laser is represented in the form 

E{z,t) = \T. £„(0exp[-/{(O„?+(^„(0)]f/„(z) + c.c. [VI-1] 



where, for example, 

in a simple a cavity laser U^[z) °^ sin k^z = sin 



n 71 



and in a ring laser t/„ (z) o« exp ( i z) = exp 



n 271 ^ 



standing wave [ VI-2a ] 



traveling wave [ VI- 2b ] 



An adaptation or interpretation of Lamb's "semiclassical" or "self-consistent" laser theory as first presented in 
Phys. Rev. 134, A1429 (1964) and refined in countless other treatments. In these lecture notes we drawn 
extensively on M. Sargent III, M. O. Scully and W. E. Lamb, Jr., Laser Physics, Addison-Wesley (1974) and 
P. Meystre and M. Sargent III, Elements of Quantum Optics, Springer- Verlag (1992). 
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With this representation for the field the induced polarization can be expressed as 

P{z,t) = ^L^Pn [t] exp[-i(a)„f+(l)„(r))] U,{z) + c.c. [ VI-3 ] 

n 

where fP„ (?) is the complex, slowly varying component of the polarization of the nth 
mode. If we take the wave equation in the form^^ 

ot c ot at 

where second term is included as a means to account for cavity losses. From Equation 
[ VI- 1 ] assuming that £ „ , fP„ , and are slowly varying functions of time^^ 

-V^E ^- = - „ exp (co„ t^K)] 1^ + C.C. 

a Z I „ a z 

^ " [VI-5a] 

= —Y,Ql exp[-/(a)„M;?+(^„)] U^ + c.c. 

n 

where Q.^ = c are the eigenfrequencies of the cold resonator eigenmodes. 

|^^^^7^ = ^L{£„-''£>„+(^„)}exp[-j(03„?+^j]t/„(z) + c.c^ [VI-5b] 



In this formulation we are assuming that V P » . 

The so called slowly-varying amplitude and phase approximation (SVAP) is used extensively in 
treating problems in laser dynamics. In the SVAP approximation it is assumed that 
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+ C.C. 



[ VI-5c ] 



3^ 3^2 = 2^ |'E„-/2 2;„(to„+^„) + '£„[-(to„+(^„) -/^„J|exp[-/(a)„?+^j]c/„(z) + c.c. II-5d] 



which reduces in the S VAP approximation to 

32 E d^E{z,t) 1 



7^ -:;L{-i2£„a)„ + £„[-a)^-2a)„(^j}exp[-/(a)„M;?+(t)„)] C/„4z) + c.c. [VI-5d'] 



and from Equation [ VI-3 ] 

a'p d'p{z,t)_i 



^ = 2^i^''>" {(0„+(^J + ^P„[-{(0„+(^„) - i(i)„J| exp[-i (a)„?+^J] f/„(z) + c.c. 
{-«n ^P„}exp[-/ (o)„?+(^„)] f/„(z) + c.c. 



[VI-6] 



If these representations are to be valid when substituted into Equation [ VI-4 ], the 
following equations must hold: 



n n 



m £ - i 2m £ -co' £ -2m 6„ co' ^„ 

n n n n n n n in n ^ n n 



[VI-7] 



We adjust \hQ fictional conductivity to account for the loss of energy or time decay of 
the given mode ~ viz. we take 0= 8„ co„/g „ where Q n is the "Q" of the mode. Equating 

real and imaginary parts of Equation [ VI-7 ] we obtain Lamb's master laser self- 
consistency equations -- viz. 
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[ VI-8b ] 


or a)„+(^„»r2„ 


2e<, ^ 


}/£„ 


[ VI-8b' ] 



Polarization of the medium ; 

The simplest application of Lamb's theory assumes that the laser medium consists of an 
ensemble of two-level atoms with a single well-defined transition ^(0^;, which is 

homogeneously broaden . Possible "pump" and "decay" processes are indicated 
schematically in the following figure: 




The model describes, in effect, a so called four-level homogeneously broaden laser 
where both levels relaxed to some presumed ground state and are "pumped" via higher 
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energy excited states, p [s,z,t',t) is defined as the density matrix operator at a place z and 

a time t which is associated with atoms excited to the state s ( where s = a or & ) at a time 
f and 'k^[z,t') as the rate at which atoms are excited to the state s (per atom, per second 

per unit volume). The macroscopic polarization is, thus, given by 



P{z,t) = Y,l df l,{z,t)(e mp ) 

sw 

= pY^j dt'\{z,t'}p,^{s,z,t\t) 



[VI-9] 



+ c.c. 



If we project this expression on to the modes of the laser cavity and compare with Equation 
[ VI-3 ], we see that 

fP„(0 = 2 ^exp [ / (a)„ t+,^it))] I dzu: (z) l^jdt'Kizj) p., {s,z,t\ t) [ Vl-lOa ] 



where ^= \ dz \ U^[z) \ is the mode normalization factor. To facilitated the integration 



of this equation, we re-write it in terms of the population matrix operator ^ {z,t) — 



VIZ. 



^it) = 2 pexp [ / (q)„ t+K{t))] ^]dz K (z) (z, t) [ Vl-lOb ] 



n 



Where p {z,t) ^Y^l dt' K{z,Op {s,zy,t) ■ [VI-U] 

S -oo 

The equation of motion of this population matrix is found by differentiating its defining 
equation ~ i.e. 
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=LK{z,t)p{s,zxt)+I,l dt'\{z,Op{s,z,t',t) . [Vl-12] 

s s ^ 

By definition p^p (5,2, =6 „j 6p 3, so that first term on the right-hand-side is replaced by 
the matrix 

X,{z,t)\ 

If we interpret the p {s,z,t',t) components of the second term on the right-hand-side of 
Equation [ VI-12 ], as the time derivatives of the pure state density matrix components 
given by Equations [ III- 14a ] and [ III- 14b ] of this set of lectures notes,^^ the component 
equations of motion for the population matrix p {z,t) become 

Paa = - Ya Paa " '^'ab f^ab + Cc] [ VI-13a ] 

Pbb = ^b - Pbb + [i '^ab Pab + c.c] [ VI-13b ] 

Pab = -('«ab+Yab)pab + '^"'l^ab [paa-Pbb] [ VI-13c ] 

Formal integration of Equation [ VI- 13c ] yields 



Including damping, Equations [ III- 16a ] and [ III- 16b ] of this set of lecture notes would have the form 

Paa = -YaPaa-[^^"''J^abPab + C-C.] 
Pbb = -YbPbb + [i^"''^'abPab + C.C.] 

Pab = 0)ab +Yab) Pab + « ab [p aa " P bb ] 
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Single Mode Operation; 

For a single-mode standing wave field in the rotating- wave approximation 

'^ab(2. =-\p'^Xt) exp [-i {«„ t'+^XtD] U. {z) [ VI-15 ] 

^ / ph-' cxp[-io)j] (z) J '^J/exp|-[/(a),,-(o„)+Y,b](?- t)- i^,{t)^'E ,{t) [paal^-O- Pbb(z>0] 
^ipff' Qxp[-i(xij] U^{z)\^ t/Texpj-[/(a),,-a)„)+YjT-/(^„(?-T)|£„(?-T) [M^'^- '^)] 

The integration can be done simply, if the changes in £ (|) and 1^^,^^ - Pj.^ are negligible 
in a time so that 

Pab(2'^) = -^ '"^^^'^ '^.{t) exp[-/(c0„ „(?))] C/„(Z) [Paa(Z'^)-M2'^)][^"Kb-«„)+Yab] ^ 

^ [VI-17] 
= -- ipr exp[-i(a)„ [paa(^>?)-M^'^)]^^Kb-«n;Yab) 

where 2^(m- viw) = [i (m- v) + w] represents the so called complex Lorentzian 
denominator. Substituting this result into Equations [ VI- 13a ] and [ VI- 13b ] we find 
the previously discussed rate equation approximation for the population components - 
- viz. 
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P'aa= ^a-YaPaa-^(Paa-Pbb) 
f^bb= ^b-YbPbb+^(Paa-Pbb) 



[ VI-18 ] 



with the rate constant is given by 



2 



^^£„^(0|f/„(z)fY.,[K-cojVYL|' 
^ k„^(r)|f/„(z)fY-;,x(a),,-a)„;Y,,) 



[ VI-19 ] 



where X(m -v; w) = [(m -v)^ + w^] represents the dimensionless Lorentzian 
function. 



For steady state ~ i.e. fj^^ = Pbb ~ ^ 



Paa-Pbb 



?^aY:-^bY;' N{z) 



Ya Yb 
lYa+Yby* 



\ + R{z)lRs 



[ VI-20 ] 



where N{z) is the unsaturated population difference and ^5 = Ya Yb (Ya+Yb) 
Combining Equations [ VI- 10b ], [ VI- 17 ] and [ VI-20 ] we obtain 



N{z)\USzt 



R{z)/Rs 



[ VI-21 ] 
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Although this result can be integrated exactly, there is more insight to be gained by 
expanding the safara^/on denominator in powers of R{z)lRs — viz. 



^ii) =- ^"'£„ 2^co,b -«n; Yab) ^]d^ M^) I UXz) f £ 



[ VI-22a ] 



To second order in R[z) /R^ we obtain 



4 " 



Ya+Yb 



V lab y 



4^ab-(0n;Yab) 



[ VI-22b ] 



where 



and 



M')-7jj-J*Mz.')|(^.W 

n 



^^l^^E/ 1 



" 2 Y,Y, 



[ VI-23a ] 



[ VI-23b ] 



is the dimensionless mode intensity. 

To complete the self-consistency loop, we combine our expression for the complex 
mode polarization. Equation [ VI-22b ], and the self-consistency conditions. Equations 
[ VI-8]. 

From Equation [ VI-8a ] 



27 Exactly 



where 



f[w) = - l-(l + w)" 



and w= T 



f \ 

Ya+Yb 



4^ab-(0n;Yab) 
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2 " ne„ 



Ya+Yt 



4^ab-«„;Yabw) 



J 



Multiplying this equation by 9 ^ ^ , we obtain 

lb lb 



I 



CO 

n 



1 ^ 

2 "he^ 



' O - 1" 3 

— Jx(co,,-co„;yJ A^^l--/„ 



^VKb-co„.;Yab)] 
^ Yab ; JJ 



[ VI-24 ] 



[ Vl-25a ] 



This nonlinear equation of motion may be written to advantage as 



[ VI-25b ] 



where the linear gain is measured by 



1 . 



2 "he^ 



V'ab J 



4^ab-(0n;Yab) 



[ Vl-26a ] 



and nonlinear saturation by 



8 ^e,, Jab 



Ya+Yt 



^l0)ab-«„;Yab) 



[ VI-26b ] 



From Equation [ VI-8b' ], we see that 
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Thus 



a)„+(|)„=n„ + -a)„^ ^ , — ^i:(co,b-co„;Yab 



" 2 

= ^„+Cy„-Pn^„ 



Yab 



Ya+Yb 



^(«ab-«„;Yab)| 



[ VI-27 ] 



where the linear mode pulling effect is measured by 



"""^ 2 ^ ITT ~1 — '^K-''-^-) 

^ "*^o lab lab 



[ VI-28a ] 



andnonlinearpMS^/ng effect by 

p _3 a)„p^ Kb-con) 

" 4 ^8„Y,b Yab 



Ya+Yb 



[4^ab-«„;Yab)] 



[ VI-28b ] 



From Equation [ VI-25a ], we find the "at-resonance"threshold condition to be: 



2QJ 2 "M^ 



ab J 



1 p^N, 



Th 



On Ky 



ab 



[ VI-29 ] 



where A^^^, is the required population inversion at threshold. 



For the general steady state condition , we have 



/ = - 



2 r Ya Yb 



- «„/P„ 



[ VI-30a ] 
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or 



4 

I.= - 

" 3 



Yab ^4^ab-»n;Yab)- A^Th/A^ 



Single- MODE, Standing- wave Laser Intensity 



[ VI-30b ] 




The time dependent intensity build-up follows from integrating Equation [ VI-25b ] 
to obtain 

«n /n(0)exp(a„r) 



a„+p„/„(0)[exp(a„f)-l] 



[ VI-29 ] 
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Transient Intensity Buildup 




2 4 6 S 10 

time 



MuLTiMODE Operation ; 

When we consider the possibility that two or more laser modes may be simultaneously 

excited or oscillating, things get quite a bit more complicated and we must refine our 

analysis. In particular, we must treat the time dependence of the population difference 
Paa ~ P^hb with greater care since variations in the population difference tend to modulate 

and couple the possible modal excitations. 

Reflection on the form of Equations [ VI- 13a ] and [ VI- 13b ] in the absence of modal 
excitation, suggests that we take 

N{z,t)^[pT {z,t)-i^ll\z,t)]= X, y-J-K r: [VI-31 ] 
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as the zeroth-order approximation to the time dependence of the population difference. 
With finite excitation of one or more modes we may write in the rotating-wave 
approximation 

'^M'^ =- ^pL2;„(0 exp[-/(a)„ U,{z) [ VI-32 ] 

^ n 

and use Equation [ VI- 14 ] to obtain the following first-order approximation to the time 
dependence of the off-diagonal component of the population matrix: 

= / n-' J j/exp[-(/a),,+Y,,)(^- o]'^a.(^,o[p:H^>^)-p:n^,^)] 

=-^/^r' A^(z,?)££„(?)exp[-/((0„?+(^„(?))]f/„(z) Jj/exp{-[/((0^^ [VI-33] 
=-Uph-' N{z,t)Y,1^Xt) exp[-/(co„ ?+^„(?))]f/„(z) 2)((0,b-(0„;Yab) 

^ n 

In carrying out this integration we have, again, assumed that changes in £ ^, and (^^^, are 
negligible in a time y~l . Using this result and Equation [ VI- 13a ] we obtain the time 

derivative of the second-order approximation to the time dependence of the population of 
the upper level ~ viz. 



r aa r aa r aa 

= ^.-Y.P..- 



2h 



N{z,i) £££ „ exp {/ [(a)„ -co J?+^„ -(t) J} Uj^z) ^03^^-(0^; y J + c.c. 



1 [ VI-34 ] 



The key point here is that the population has pulsations at the intermode beat 
frequencies! Integrating we find 
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2 

N{z,t)YJl'£.'£r. exp{ i[K-(oJ t+<^,-<^Ji}u^{z)ul{z) 2>(co„-co„; Ya) ^D^b -«m;Yab) + c.c. [ VI-35 ] 

n m 

and since P'bb' =-pl^' with =>Yj, we obtain the second-order approximation to the 
population difference as 

2 

N{z,t)Y.'L'£.'£r. exp{/[(to„-a)J t+^,,-^Ji}u^{z) U*{z) 

n m [VI-36] 

x[®(a)„-to„;Yj + ®(co„-a)„;Y,)][2)(a),,-a)„;Yab) + ®(o>„-a),,;Y,J] 

Substituting this expression into the formal integral of Equation [ VI- 14 ], we find, 
directly, the third-order approximation to the off-diagonal component of the population 
matrix 

2 

N{z,t) IH,^{E^'E.Uiz) K{z) U,{z) 

1 m n 

xexp{-i[(03 -(o„+(oJr+^^-^^+^J}^D((o,b-cOl+03„-(o„;Yab) [VI-37] 

x[2)(a)„-a)„;Ya) + 2t«m-«n;Yb)][^Z^«ab-»n; Yab) + ^^{»m-«ab;Yab)] 

Using Equation [ VI- 10b ] and invoking Equations [ VI-8 ] the general self-consistency 
conditions may be written ^8 



r aa 



in 



K aa K bb 



in 



Pab 



in 



28 It is important to note, that use of Equation [ VI- 10b ] requires evaluation of the following integral : 

L I L 

Numn = \dz N(z) U;{z) mz) Ujz) U„iz) ndz \U,{zt 





R. Victor Jones, March 16, 2000 76 



The Interaction of Radiation and Matter: Semiclassical Theory Page 77 

£\ = «k £ k - LLLe i-E „ Im{^kin,„ exp [ / [ VI-38 ] 

1 m n 

%+<^k=^k+<?k-LLL£k'£i£m£„ Re{i^kim„ exp [ / %,^,]] [ VI-39 ] 

1 m n 



with the coefficients summarized in the following table for a standing wave 
configuration: 



Coefficients 


Significance 


«k= 4^ab-»k; Yab) -2 (^JQ^ 


Linear net gain 


<^k= [Kb-«k)/Yab]4«ab-«k;Yab) ^ 


Linear mode pulling 


Pk=^'Kb-»k;Yab) ^3 


Self saturation 


^kl.„ = i^\ r ^1 {l+ [A^.(_) + A^2(»-l) ]/^} ^(«ab -«„; Yab) 


General saturation term 


X [^2^«ab -»„; Yab) + 2^«m -«ab; Yab)] [^X^m -«„; Ya) + Yb)] 





Thus, the coefficients in Equations [ VI-38] and [ VI-39 ] are sensitive to the particulars of the laser mode 
configuration. Since the modal excitations cannot change rapidly , the only integrals of importance are those 
forwhichCO^ -CO; -l-CO^ -a)„ =0 
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t3^ +t3^ = 7 i 

"kkmm^^kmmk * ^ 



Ya+Yb 



pim f {2+ A^,(„_,)/ A^} ^D(a),, -co,; y,,) 



L YaYb 

X{4^ab -«k; Yab) 4^ab Yab) 



+ ■ 



Y, 



ab 



YaYb 



.Ya+Yb, 



[^^«ab k; Yab) + ^^«m -» ab : Yab)] 

X [ 2>(0>m - « k; Ya) + 2^0>m - » k; Yb)] } 



= [(2^'YaYb)/lrf] Im(l^tkmm+1^kmmk) 

= 1^3 {2+ A^2(m-.)/ A^} [4«ab -«k;Yab) ^^ab Yab) 



+ 



2 [Yab] 



YaYb 
Ya+Yb 



Re|^)(a),b k; Yab) [^Wab k; Yab) + ^(Om ab; Yab) 

X [!Z)(a)„ -CO,; Ya) + 2)(co„ -co ,; yj]} 



Stationary coefficient 



Cross-saturation 



%m„ =K-«l+«m-«„) ^+^k-^l+^m-^„ 
Pm = (»ab m)/ Yab ^ ' («ab Yab) ^3 

T^km = [(2^^ Ya Yb)/|pr] +1^ 

kmmk ] 

= 1^3 {2+ A^2(n,-k)/ A^} [(Wab - Wk)/Yab ^ («ab -«k; Yab) ^ (^ab Yab) 



-^[Yabf Im|^D(co,,-co,;Y,,)[25(co,,-co,;Y,,) + ^D(co^-co,,;Y,,)] 

. Ya + Y b 



Relative phase angle 
Self mode pushing 



Cross pushing 
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First-order factor 

and 

Third-order factor 


A^2m= , 1 dzN{z) COS 


\ 1 

z 

1 J 


Spatial factor 


CD^u- v;w) = [i (m- v) + w] 




Lorentzian denominator 


L{u— v; w) = w^ [( ~ + ~] 


-1 


Dimensionless Lorentzian 



Two-Mode Operation; For two modes the amplitude detennining equations ~ i.e. 
Equation [ VI-38 ] - reduce to 

[VI-40] 

2 — 2 1 ^ 2 P 2 "^2 ^ 21 "^ll 



and the frequency determining equations ~ i.e. Equation [ VI-39 ] ~ reduce to 

[ VI-41 ] 



^2+^2 =^2+^2~p2 '^21 A 



W 



Multiplying the amplitude equations by V aY^) obtain the equations of motion 

for the dimensionless intensities 

i, = 2I,[a-^J-%,,I,] 
1^ = 2 I^la^-^^ I2-Q21 A] 
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Stability of possible steady state solutions ; 

Stability criterion for stationary solutions: If /i=//*'+8 and /2=/2^'+e2, ^ and 
82 ^ as t^oo. 

Solution 1: = and 7^'' = a2/p2 so that 

£1 = %][«!- Pi 61-612 (fl2/p2+e2)] 
E2 = 2[a2/p2+e2][«2-p2(«2/p2+e2)-02i ej 



[ VI-43a ] 



or 



82 =-2 (a2/|32) [|32 £2 + 821 eJ+0(8^) 



[ VI-43b ] 



Stability requires that a^= ai-9i2 «2/P2 remain negative. If a{<0 but a^>0 one says 
that I2 inhibits the oscillation of /j (mode inhibition). If becomes large to overcome 
the inhibiting effect of 4 , the solution becomes unstable and /j builds up. 

Solution 2: 9^ and 7^'' 9^ so that 



which has the solution 



a, = 13, 71" + 9,2//' 

«2 = P2 4'' +021 A*'^ 



^W_ (^l/Pl)-(e2l/Pl)(«2/P2) _ ai7Pl 

' 1-C i-c 

^M_ (^2/p2)-(ei2/p2)(^/Pl) _ ^27P2 

' 1-C 1-C 



[ VI-44 ] 



[ VI-45 ] 
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0^, 

where the coupling constant C= . Substituting this solution into Equation 

Pi P2 



[ VI-42 ] we see that 



In matrix forni 



d_ 

dt 



^1 =-2 [pi 81+612 £2] + 0{8') 

£2 = - 2 / [p2 £2 + e 21 81] + 0(8') 



-2 r «; a;(ei2/i3i)i 

l-c[a^e2i/p2) a'^ J 



[ VI-46a ] 



[ VI-46b ] 



The solutions will be stable if the eigenvalues of this equation are both negative. 

^ Det(5^^ -Xl)=0 [ VI-47] 







= X 












£2. 




£2. 



^i + aj 



(■.)■ 



1-C 



1-C 



-4 



^1 ^ 2 

vl^y 



[ VI-48 ] 



Case 1: 



Case 2: 



Weak coupling where C<1 and ai'>0 and a2>0: Stable since 
both eigenvalues negative 

Strong coupling C> 1 and a[<0 and a2<0: Unstable since one 
eigenvalue is positive 
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